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We investigate mechanisms of the typically observed recoverable prevalence in epidemic spreading. 
Assuming the time-independent connectivity correlations, we analyze the dynamics of spreading on 
linearly growing scale-free (SF) networks, and derive the extinction condition related to the rates 
of network growth, infection, and immunization of viruses. The behavior is consistent with the 
previous results for SF networks by a mean-field approximation without connectivity correlations. 
In particular, it is suggested that the growing must be stopped to prevent the spreading of infection. 
This insight helps to understand the spreading phenomena on communication or social networks. 
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In the world-wide communication and transportation, 
we arc able to work on rapid consideration or decision 
for scientific, technological, social, and business trends, 
while we are also frighten by threat to dangerous viruses 
for computer systems and our body itself. The needs 
become larger for understanding the mechanism of epi- 
demic spreading and for finding effective strategies of 
the prevention against huge damage in economy or life 
by computer viruses via e-mails [10] [21], human im- 
munodeficiency virus (HIV) [12], severe acute respira- 
tory syndrome (SARS) [Ki], and so on. For the pur- 
poses, the study of infection disease on network models 
is useful. In particular, we should remark that many 
real networks, whose node and link are corresponding 
to individual (computer or person) and the (communi- 
cational, pathological, or sexual) contact between them, 
have scale-free (SF) structures based on natural rules of 
network growth with new nodes and preferential attach- 
ment for linking as rich-get-richer phenomenon [2] [4], 
in which the degree distribution exhibits a power-law as 
P{k) ^ k '', 2 < 7 < 3, for the probability of k connec- 
tions in address-books of e-mails [8] [9] or sexual partners 
[12]. ^ 

It is very possible that there exist a common mecha- 
nism for epidemic spreading in realistic SF networks, and 
that the extremely heterogeneity consisted of many nor- 
mal nodes with less connectivities and a few hubs with 
much connectivities is crucial for the superspreading 
and the persistence of viruses in long-time [21]. Indeed, 
the important properties of SF networks are the robust 
and vulnerable connectivities against random and tar- 
geted attacks, respectively [1]. The property is applica- 
ble for the targeted immunization for hubs [7] [19] to pre- 
vent the spreading of infection. Thus, the universal and 
robust structure in growing complex networks attracts 
physicists, and surprising results have been recently ob- 
tained. 

In striking contrast with the usual models for epidemic 
spreading [.'>], it has been shown [ I T] that a susceptible- 
infected-susceptible (SIS) model on heterogeneous SF 
networks has no epidemic threshold at a large network 



size nearly infinity; infection can be proliferated, what- 
ever small infection rate they have. This result is closely 
related to the bond percolation problem. The epidemic 
threshold in a susceptible-infected-recovered (SIR) model 
coincides with the threshold for the bond percolation 
problem on networks [15]. On the other hand, in a closed 
system of the conventional susceptible-hidden-infected- 
recovered (SHIR) or SIR model, the number of infection 
is initially increased and saturated, finally converged to 
zero as the extinction. The pattern may be different in an 
open system, in fact, oscillations have been described by 
a deterministic Kermaek-McKendrik model [20]. How- 
ever a constant population (equal rates of the birth and 
the death) or territorial competition has been mainly dis- 
cussed in the model, the growth of computer network 
or world-wide human contacts is obviously more rapid, 
and the communications of mailing or the transporta- 
tions are not competitive. Even in a new report for mod- 
eling the SARS epidemic [13], the mechanism of widely 
rapid spreading is not clear enough, though further treat- 
ments for the heterogeneity in transmission through su- 
perspreaders has been pointed out and the effectiveness 
of quarantine has been discussed. The strategic quaran- 
tine in both the model and experiences in Hong Kon is 
potentially related to prevent the typically observed re- 
coverable prevalence [ ' )] [21] in an open system as grow- 
ing with new individuals. Fig. 1 shows an example of 
typically observed pattern in the prevalence of computer 
viruses. 

Thus, we consider a growing system on heterogeneous 
SF networks for e-mails or human contacts, and inves- 
tigate the mechanisms of recoverable prevalence and ex- 
tinction of infection. In the previous paper [9], we have 
derived the conditions of extinction for the spreading in 
deterministic SIR models on SF networks by a mean-field 
approximation without the connectivity correlations. Al- 
though the connectivity correlations are not found in all 
growing network models and real systems [18], they are 
at least quantitatively significant for the spreading on 
heterogeneous SF networks [5] [G] [11]. In this paper, we 
extend the previous results to correlated cases between 



2 




2 3 4 5 
2003 

source: i/^'lc. 1rer-clnicrci.com 



FIG. 1: Typically observed pattern as recoverable prevalence 
of the computer virus NIMDA extracted from [21]. 



the degree of connections, and discuss the mechanisms of 
recoverable prevalence and extinction on linearly growing 
SF networks. 

Heterogeneous SIR model.— On heterogeneous SF net- 
works, we consider a deterministic SIR model [14] with 
macroscopic equations for the number of states. The 
state transition is from susceptible, infected, to recovered 
or immune, whose numbers for each connectivity k are 
denoted by Skit) > 0, /^(t) > 0, and Rk{t) > 0, respec- 
tively. Microscopically, susceptible nodes stochastically 
become the infected by contacts with infected nodes, 
which are recovered by the detection or immunization. 
Assuming that infection sources exist in an initial small 
network and that both network growth and the spread 
of viruses are simultaneously progressed, we introduce a 
linear kernel [11] as Nk{t) ^ x t, Nk{t) = Sk{t) + 

Ik{t)+Rk[t), the growth rate Ofc Ak"", A>0,v>2. 
Note that the total N{t) = f^k^kit) ^ (Efcafc) x ^ 
means a linear growth of network size. Since the max- 
imum degree increases as progressing the time and ap- 
proaches to infinity, it has a nearly constant growth rate 
Sfelm «fc ~ Im Ak-^dk = ^^fzT- for large t. As shown 



in [11], the introduction of linear kernel is not contradic- 
tion with the preferential (linear) attachment [2] [4] . 

We first consider a simple case with only the detection 
of viruses. The time evolutions of S'fc > and Ik > for 
each class of the connectivity k arc given by 



dSkjt) 

dt 
dlkjt) 

dt 



-bkSk{t)Qk{t) + ak, (1) 
-5oIkit)+bkSk{t)ek{t), (2) 



where b and 5a denote the infection and detection rates 
between and 1, the shadow variable Rk{t) is implic- 
itly defined by = Sohit). The factor Qkit) '^^^ 
J2i ^Ii{t) represents the expectation of infection by 
contacts from degree classes {1} to the degree class k. 



We assume the correlation between degrees Uki ^ fiik, 

n-k X); "-fcij ^rid Nki ~ nki x t which is defined by the 
number of nodes with degree k attached to a node with 
degree I at the time t. We can easily check = Uk and 
the solution Nk{t) ~ Nk{0) + ak x t ^ ak x t because of 
the initial small size (almost all Nk{Q) = except for a 
few small k). The network is keeping the power-law de- 
gree distribution P{k) A:"'' in growing. Although the 
growing in our model is not exactly same in Krapivsky- 
Render's GN model [11], the growth rate is close to 
their rij. = ,,,,, , ^f,,, , ^ 4fc~'^ in large t. 
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Fig. 2 (a) shows the nuUclines of 
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and the vector filed for Eqs. (1)(2). 
equilibrium point: It 
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bk'-^i^+i ^ ukii^'-' ' state space, the viruses are not 

extinct by only the detection in even /* = (Z ^ fc), 
unless the growing is stopped as Ofc = 0. 

Effect of immunization.— Next, we study the effect of 
random and hub immunization. With the randomly im- 
mune rate < Sr < 1, the time evolutions are given by 
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-bkSkit)ekit)+ak~SrSkit), (3) 
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where the shadow variable Rk{t) is also defined by 

^=SaIk{t)+Sr{Sk{t)+Ik{t)). 

On the section by assuming 31^ > ior I ^ k, the 
nuUclines are 
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for Eqs. (3) (4). At the intersection of nuUclines, from 
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When the right hand side of (5) is denoted by /(O^), for 
3Qk > 0, ^^^^g^|efc=o > 1 is necessary. The condition is 
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FIG. 2: NuUclines and the vector fields near an equilibrium 
point, (a) damping oscillation on a section /;*: const, for all 
I ^ k, (b) nearly extinction on a section /,* « 0, (c) divergence 
from a saddle, and (d) convergence to the extinction 7^ ^ 
on a section /* =0. 



This means the growth and infection rates should be 
larger than the immune and detection rates according 
to the connectivity correlation of degree k. 

On the other hand, we assume /;* = for I ^ k to dis- 
cuss the extinction. The necessary condition of extinction 



is given by that the point (0, j^) on the hyperbolic nuU- 
is below the hne Sk = Sk =^ 



cline of ^ 
const, of ^ = 0. From the condition 



Ok ^ (^0 + Sr)nk 
5r khukk 



we obtain 



5r > -So + \ /5^ +4fcafcb x n kk/nk 
~ -So + ^yS^ + Akb X Ukk, 



(7) 
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for Uk ttk asymptotically. From Sr < 1, we further 
derive 



A < min 

k>rn 



(1 + 2So)nkm'' 
bukk 



as a limitation of growth rate or a weak correlation 

l + 2So 



nkk < 



Akb 



for the extinction. In this extinction case, there exists a 
stable equilibrium point, otherwise a saddle and a stable 
equilibrium point as shown in Figs. 2 (c)(d). Note that 
Sr > is important for the extinction. When no growing 
Qk = (as a closed system) is applied to the above results, 
the condition (6) for damping oscillation is not satisfied, 
while the extinction condition (8) is satisfied. In other 
words, it is suggested that a growing network causes the 
recoverable prevalence of infection. 

By replacing Ofc = Ak~'^ with the Krapivsky-Render's 
Uk [11] and using their Ukk = k^kiiyH2k+t)%k-i) ^ ^ ^ 2, 
we can check that Ukk < ^l^kb" satisfied for any < 
6, (5o < 1 . The reason why viruses tend to be extinct 
in their GN model by even random immunization is that 
the infected parts are disconnected by recovered nodes on 
the tree structure, while our model generally has other 
bypass links than the links of tree in the difference of 
the degree distributions akt and their rikt for small fc, 
therefore infection can be spreading. 

Moreover, for the hub immunization [7], Sr is replaced 
by < Shk'^ < k'^ , r > 0. Then, the necessary condition 
of extinction in (8) is relaxed to larger 
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However this is under the assumption of 
general case for a section /;* ^ is similar to Fig 
(a), then viruses in the class of degree k are not extinct 
unless the growing is stopped as = 0, otherwise they 
approaches to the extinction as shown from Fig. 2 (b) 
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to (d). In fact, at the point (. . . , I^, Sk, ...,/;, S*;, .. .) 
(. . . , 0, 1^, . . . , 0, f^, . . .), the eigenvalues of Jacobian 
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are A = < and A 
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< from the condition (7), 



+ fe^^^ and Cki 



^^^HM. ^ if (;^^ and with respect to the cross correlations 



are negligibly weak. Another extinction route is given by 
6 0, in which case the condition (8) is satisfied for any 
A and Ukk ■ It is consistent with the result of no epidemic 
threshold [17]. 

These dynamics are qualitatively same to the previ- 
ous results by a mean-field approximation without the 
connectivity correlations [9] . In summary, we have inves- 
tigate mechanisms of recoverable prevalence and extinc- 
tion of viruses in SIR models with the connectivity corre- 
lations on linearly growing SF networks, and derived the 
necessary condition of extinction related to the growth, 
infection, detection, and immune rates. The results sug- 
gest that we must stop the growing of network or close 
it by quarantine to prevent the spreading of infection. A 
more general case in linearly growing SF networks on no 
assumption of correlation Nki ~ Uki x < is further study. 
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